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We study the quantum-mechanical interpretation of mod- 
els with non-Hermitian Hamiltonians and real spectra. We 
set up a general framework for the analysis of such systems in 
terms of Hermitian Hamiltonians defined in the usual Hilbert 
space Z/2(— oo, cxd). Special emphasis is put on the correct 
definition of the algebra of physical observables. Within this 
scheme we consider various examples, including the model re- 
cently introduced by Cannata et al. and the model of Hatano 
and Nelson. 
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I. INTRODUCTION 

Recently, models with non-Hermitian Hamiltonians 
(non-Hermiticity meant here in the sense of the space 
L2(— 00,00)) have attracted a lot of interest. One of 
the earliest attempts to use such systems is the work 
of Hatano and Nelson Q where it was suggested that a 
delocalization transition in superconductors can be de- 
scribed by a non-Hermitian Hamiltonian. Next, Bender 
and coUaborators in a series of papers |^ investigated 
some non-Hermitian, VT symmetric Hamihonians, ar- 
guing that they have real eigenvalues. This has triggered 
a lot of activity in this field ||^-[ll|| . 

Despite this, some fundamental issues concerning the 
quantum-mechanical interpretation have only started to 
be addressed p^-p^. In the models treated by Bender 
et al. it is, for example, necessary to extend the defini- 
tion of position-space wave functions to complex values 
of the coordinate . This means that the wave func- 
tions are not elements of the Hilbert space L2(— 00,00), 
so that the notion of non-Hermiticity in the sense of the 
space L2 does not seem to be useful here. In our opinion, 
the physical meaning of these wave functions deserves 
additional study. 



In this paper we will address these questions. In Sec- 
tion U we start by constructing a Hilbert space 7i that 
contains superpositions of eigenfunctions with real eigen- 
values of a Hamiltonian H. Here the choice of a scalar 
product that provides the link to a probabilistic interpre- 
tation is quite arbitrary. Motivated by the approach sug- 
gested by Bender and collaborators [||, we consider, e. g., 
scalar products that are defined along complex paths. 

We remove the arbitrariness for the choice of the scalar 
product by demanding that the Hamiltonian is to be 
interpreted as the generator of the time evolution for 
a closed system, and therefore has to be Hermitian in 
the space Ti.. This still leaves some arbitrariness for the 
structure of Ti, but all allowed theories have a consis- 
tent quantum-mechanical interpretation. Next we con- 
sider the algebra of physical observables in the theory. 
We set up a canonical formulation in which the Hamil- 
tonian is only a function of two Hermitian operators x'^ 
and p'^ fulfilling canonical commutation relations. Such 
operators may under some conditions have an interpreta- 
tion as a position-space variable and the generator of its 
translations, resp. One can then employ the uniqueness 
theorem by von Neumann ||l^ to map the space TL to the 
space L2 ■ In this way, the physical interpretation is fixed 
uniquely. 

The whole constructi on is applied to two exactly solv- 
able examples in Section HI , The first example is a simple 



non-Hermitian model based on the well-known, Hermi- 
tian one-dimensional Coulomb problem on the real half 
axis. Here we show that the canonical formulation does 
not contain any new information. The second example 
is based on the model introduced by Cannata et al. IT^l 
in which the potential is given by V{x) = e 
we also construct a canonical formulation H 
with Hermitian H, but then we show that this formu- 
lation does not allow the interpretation of and as 
position and momentum operators, resp., so that the in- 
terpretation of the model remains unclear. We end this 
Section with a discussion of the phenomcnologically im- 
portant model of Hatano and Nelson 0] in the spirit of 
our approach. 
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Our concluding remarks are contained in Section IV 
Some technical details that are used in Section III are ex- 
plained in Appendix In Appendix |^ we briefiy com- 
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ment on the connection between real spectra and VT 
symmetric Hamiltonians. 



II. GENERAL FRAMEWORK 

Suppose we have a non-Hermitian (in the usual sense 
of the space Li{—oo^ oo)) Hamiltonian H and some eigen- 
functions '\\}n{x) with real eigenvalues £'„, 



En real 



(2.1) 



Here the operators x and p act in the usual way as mul- 
tiplication by X and — i times differentiation, resp., on 
ipn (x) , and the eigenfunctions are calculated as solutions 
of the corresponding differential equation. They are not 
required to be normalizablc with respect to the norm 
11^11^2 — y/{ip,ip)L2 of the Hilbert space L2(— 00,00). 
(We use the notation (V', ■(/'')i2 = dxi/'*(a;)V''(a;) for 
the scalar product of L2-) 

Our aim is to set up a formulation in which we can 
interpret the superpositions of the eigenfunctions ipn 
quantuni-mechanically. For simplicity, we will do this 
only for the case of a discrete, infinite, non-degenerate 
spectrum of H. 

To start, we define the vector space V = sp8Ln{ip„ , n — 
0, . . .} of the finite superpositions of the eigenfunctions. 
Since eigenfunctions corresponding to different eigenval- 
ues are always linearly independent, the dimension of 
V is infinite. Next we define a scalar product (., .)v 
on this space. Besides requiring the usual properties 

((^,V)t„ II VII?. = (V',V')v > for V and Uh = 
for ^ — 0), we leave this scalar product arbitrary at this 
point. The scalar product turns V into a separable Eu- 
clidean space, and we can use standard theorems JT^ to 
complete this space, thereby defining a separable Hilbert 
space 7i. As a result, we have V d Ti. and 



(^, V^')^ = (V>, V/)v for all V, V^' e V 



(2.2) 



(Recall that in addition to the finite superpositions of the 
Ipn, n = 0, . . ., that make up V, the Hilbert space H also 
contains all limits ip = lim„_+oo/n, fn G V, of Cauchy 
sequences of vectors of V.) 

The Hilbert space H is the natural choice as the space 
of states for the system described by the Hamiltonian H, 
because it allows a consistent probabilistic interpretation 
of the model in terms of scalar products of states. 

It is well-known ||l^ that all infinite-dimensional, sep- 
arable Hilbert spaces are unitarily equivalent to the 
Hilbert space i2(— oo,oo). This unitary equivalence 
means that an isomorphism T : 7i — > L2 exists that re- 
spects the scalar products in both spaces: 



from a unitary operator U that is defined to be an au- 
tomorphism, e. g. U : L2 L2, that respects the norm, 
and therefore fulfills = U^^.) Given an operator A-j-c 
that is defined in Ti, there is a corresponding operator 



that is defined in L2 and fulfills 

{i;,An^')n^{Ti;,AL,T,p')L, 



(2.4) 



(2.5) 



for all tp Ti, and ip' in the domain of definition of A-j-c- 
If one denotes the Hermitian adjoint of an operator An 
with respect to the scalar product in Ti by A^^, 



finds 



{i;,Ani^')n = {Al,i;,i^')n 



a{^ = taI^t-'^ 



(2.6) 



(2.7) 



so that, for example, A^^ is L2 Hermitian if A-j-c is Ti. 
Hermitian, and vi ce v ersa. 



The equations (2^) - (2/7) imply that theories defined 
in Ti. and L2 are physically indistinguishable. 



A. Hilbert spaces with Hermitian Hamiltonian 

We have, up to this point, not specified the details of 
the scalar product in Ti.; the construction outlined above 
works with every scalar product that can be defined on 
the space V. But obviously not only the value of ma- 
trix elements, but the whole structure of the space Ti. 
depends on this choice. And although any choice allows 
a consistent probabilistic interpretation, we will here con- 
sider only Hilbert spaces in which the Hamiltonian is 
Hermitian. T he reason is that we consider the system 
described by (2T) as a closed system. (This appears to 
be an implicit assumption in most of the models treated 
recently.) According to a theorem by Wigner ||l^ the 
unitarity of the time evolution for such systems is a con- 
sequence of fundamental properties of quantum theories 
as, for example, the linearity of the time-evolution oper- 
ator and the fact that physical states are described by 
rays {Xip; A G C, ■0 G Ti} rather than vectors tp. If we 
want to keep these properties, and if the Hamiltonian H 
is to be interpreted as the generator of time evolution, H 
has to be Hermitian in the underlying Hilbert space H. 

In our context, it is the reality of the spectrum that 
makes the construction of Hilbert spaces with Hermitian 
H possible. Define, for example. 



{'(pn,1pm)H = Sn 



(2.8) 



{il},iP')-H = {TiP,Ti^')L,_ for all e H . (2.3) It is easy to verify that (^ 



The transformation T is called a unitary transformation 
from Ti. onto L2 . (This notion has to be distinguished 



is a well-defined scalar prod- 
uct in Ti. (Here the linear independence of the eigenstates 
Ipn is important.) We have (from now on we denote H 
hyHn) 



2 



n 



E„Sr, 



H 



(2.9) 



and since the i/'n span the whole space, we can conclude 



and 



(2.10) 



where we have used (|2.7| ). Thus, the choice (|2.8| ) leads to 
Hermitian Hamiltonians and H^^ their respective 
spaces. 

Note that there is no problem to explicitly construct 
a transformation T : Ti ^ L2 and the Hermitian oper- 
ator Hl^ . It is enough to take an arbitrary complete, 
orthonormalized set of L2 functions, e. g. the eigenstates 
^Pn{x) of the harmonic oscillator, and define T in (2.3) as 
the linear transformation that fulfills T^/'n = Vn- Then 
Hl^ — TH'n{x,p)T^^ has the required properties. The 
crucial point is that the transformed operators XL2 and 
will have the same complicated properties concern- 
ing their Hermitian adjoints as x and p and will therefore 
in general be void of any physical significance. One has 
to find a transformation T that leads to a Hamiltonian 
Hl^ — H{xl2iPL2) thai allows a clear physical interpre- 
tation. 

For achieving this, one may employ the uniqueness 
theorem by von Neumann It states that all irre- 

ducible representations of two self-adjoint operators x"^, 
p'^ that are defined in a separable Hilbert space and fulfill 
canonical commutation relations arc unitarily equivalent. 
This means that given two such operators in a separable 
Hilbert space 7i, a unitary transformation T : 7i — s- L2 



exists such that = 
canonical L2 operators 



Tx^yT-^ 



and p\^ are the two 



{x'iip){x) = X(p{x) , {p'iw){x) 



.d^J 

dx 



(2.11) 



In this way the uniqueness theorem allows one to find 
a unique (modulo unitary equivalence) interpretation for 
quantum theories that are defined on separable Hilbert 
spaces other than L2- 

Thus our aim is to find a set of canonical operators 
(a;^,p^) that are self-adjoint in 7i, 



iPn)^=Pn 



and fulfill 



(2.12) 



(2.13) 



If we are able to express H-h as a function of these two op- 
erators, Ht-i = H {x^, p^) , we can immediately transform 
the model into the space L2 with H^^ = H{x'i^,p'i^). 
Here the physical meaning is clear, because one can in- 
terpret x'j^^ as the position-space observable and as 
the momentum (generator of x'j^^ translations). 

We emphasize that this canonical formulation H^^ = 
H{x\^ ,P%^) is not just an alternative formulation for the 
original problem. In our opinion it is the only formula- 
tion which has a clear physical meaning. Thus, if such 



a formulation turns o ut t o be impossible, we argue that 
the original problem (2J^) is inconsistent. On the other 
hand, if a canonical formulation is possible, then the non- 
Hermiticity of (2.1) with respect to the space L2 is only 
a superficial one, due to a choice of variables that is in- 
adequate for the given physical problem. 

A necessary condition for the self-adjointness of the op- 
erators a;'^ and p*^ is their Hermiticity. In physical consid- 
erations both notions are usually identified. We shall do 
the same here, except in cases in which the difference be- 
comes essential for the physical conclusions. These cases 
will be clearly indicated. 

Of course, definition ( |2.§| ) is not the only scalar 
product giving a Hermitian H. The condition that a 
scalar product (., .)y in the space V leads to a Her- 
mitian Hamiltonian H follows from the requirement 



)v 



(i?V'n,V'm)v all n,m, which gives 



{En — i?m)(V'nj V'»n)v ~ ^- Thus all scalar products for 
which the ipn are mutually orthogonal lead to a Hermi- 
tian H. Therefore the most general ansatz for a scalar 
product with Hermitian H is 



(V'n,'0m)v = lnS„ 



(2.14) 



where the 7„ are some positive constants. With respect 
to this scalar product, vectors ■(/'„ = ipn/cn with |c„p = 
7„ are orthonormal. Let's define the Euclidean space 
£ to be the space V equipped with the scalar product 
, and let £ be the Euclidean space build from V and 
Eventually, £ and £ would be completed to give 

Every linear operator 



(2 



m3- 

the Hilbert spaces H and 7i, resp. 



AllJn 



(2.15) 



has a meaning in both spaces £ , £ ^ but their properties, 
for example concerning their Hermitian adjoints, will in 
general be different in £ and £. On the other hand, the 
Euclidean spaces £ and £ are connected by the unitary 
transformation T : £ — > £, 

?/;„ . . 

= VAi = — , (V'n)V'm)v = (rV'n,TV'm)v , 

(2.16) 

and therefore A as an operator £ — > £ can be mapped to 



A = TAT^' ■.£^£ 



(2.17) 



These operators, A : £ £ and A : £ £, sh are the 
same properties. The transformation T in ( 2.16| ) is only 
a function of iJ, so that H = H . But in general, unless 
A commutes with T, one has A^ A. 

This means that to a certain extent one can adjust 
the properties of an operator A : V ^ V hy choos- 
ing an appropriate scalar product. This will turn out 
to be helpful for t he co nstruction of sets of canonical 
operators (Section IIIB). If, for example, an operator 
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F{H)-'^AF{H) :£ ^£is Hermitian in then A is Her- 
mitian in the space £ which is defined by = F[H)ipn 

and {■!/'„, ?^m)v = ^nm- 

B. Scalar products along complex paths 

In the work of Bender and collaborators ||] the eigen- 
functions ipn^z) are defined with z describing a curve in 
the complex plane which is chosen in such a way that the 
eigenfunctions are decreasing asymptotically. Let a pa- 
rameterization for this curve be 2; = z(s), —00 < s < 00. 
We may then define a scalar product by 



/ ds[^(z(s))]*^'(z(s)) 



(2.18) 



The question is whether such scalar products can lead to 
a Hermitian Hamiltonian. 

According to our general discussion, is an element 
of Ti. and ip{s) = Tjj{z(s)) € L2, i. e. the transformation 
T : H ^ L2 acts as 



and 



(V,V'')H = (rV',Wi 



(2.19) 



(2.20) 



Consider two sets of operators: a; 7^ and pu act in Ti in 
the usual way: 

{xn^){z) = z^Piz) , (p„V)(^) = -i^ (2.21) 

and can be transformed into xl^ and PL2 which act in a 
complicated way. The second set consists of the canonical 
operators x^^ and p^^ that act in L2 like 

{xl^^){s) = s^{s) , (pc^^)(,s) = -i^ , (2.22) 

are Hermitian and can be transformed into a;^ and p^ . 
With these operators, one finds 

(2.23) 
(2.24) 

= if,z{xlJ'p')L2 = {'f,XL2'p')L2 , (2.25) 

i. e. XL2 = Tx-hT^^ — z(x^^), which is the path z{s) with 
the real parameter s replaced by the Hermitian operator 
cc^^. This is of course also true for the operators in Ti: 



{ip,xnfp')n = / dsip*{s)z{s)ip'{s) 

ds^*is){zixlW)is) 



xu = z{x''^) 



(2.26) 



The operator pq-i can be analysed analogously: Using 
Atl^/dz = {dip / As) / {dz / As) , one obtains 



(V',PwV'')« = -i / ds[^(z(s))] 



dV^' 



As >f*{s) 

" Az/dxt^ 
and consequently 

PH = 



z—z(s) 

Pl2^']{s) 



Az{xlJ/Axl 
Pl2'p')l2 = {'P,PL2'p')l 



dz(a;^)/dx^ 



Ph 



(2.27) 



(2.28) 



Hence definition (2.18) leads to simple relations connect- 
ing the original operators {x-h , p-u) with the canonical 
ones (a;^,p^). (But note that ( |2.18 ) explicitly depends 
on the parameterization chosen. This leaves a certain 
freedom for the definition of the canonical operators.) 

The Hermiticity of the canonical operators immedi- 
ately leads to 



(2.29) 



and with Axu / Ax'j_j)'H = p^dx^/dx^ and [pn, /(a^w)] = 
— id//da;-^ one obtains 



with 

v{xl^) 



pIl = v{x''y)[iw{xl^) + pn] 

Axt-c/Ax^ 
dx^/da;^ 



A^x\, 



(dxw/da:^)(dxt^/dx^) (dx^)' 



(2.30) 

(2.31) 
(2.32) 



One can check that (p^)"'' = Ph\ this gives the unitarity 
of v. 



v{x'^)'^ — v{x'^) 



and an additional condition: 
Av 



Ax 



H 



Ax-H 

Ax: 



(2.33) 



(2.34) 



H 



These equations are identically fulfilled for (2.31) and 

Now one can ask whether a curve or parameterization 
exists that makes a given Hamiltonian Hermitian. Con- 
sider a Hamiltonian of the form 



H{xn,Pn) =Ph + Vi^n) 
nitia 

operator (p^)^ is given by 



In its Hermitian conjugate = {pj. 



ip\i? 



AxH 



+ 1 



2v w — V- 



(vw)'^ 
Av 



(2.35) 

{V{xn))^ the 

(2.36) 



Axn 



PH + v^p'^ 



4 



Since {V{x-h))^ does not contain p-u, a necessary con- 
dition for the H ermiticity of H is that the term Hnear 
i n p'H in ( ^.36[ ) vanishes. This, together with relation 
( 2.34 ) gives 3vw — w^. Taking the Hcrmitian adjoint 
of this equation and using the unitarity of v then shows 
that necessarily w = 0. Since w contains d^a;^/(da;^)^, 
which is the operator version of d^z*/ds^, one sees that 
only straight lines, 

z{s) = a + bs , a, 6 G C , 



are allowed. Then u is a constant phase, v{xn) 



(2.37) 

b/b* = 



ari 



g*- and (py2 



arg bp2^ ^ Hcncc IS Only Hcrmi- 
tian if the parameter b satisfies 

77 



arg 6 = n- 



Under this condition, one has 

Hi=pl^ + {V{xn))^ 



(2.38) 



(2.39) 



and it may happen that a and |6| can be chosen in such a 
way that also the potential becomes Hermitian. But the 
fact that only straight lines in the comple x plan e can lead 
to Hermitian Hamiltonians of the form ( 2.35 ) seems to 
restrict the usefulness of scalar products along complex 
paths severely. 

In the simple case of straight lines one can expl icitly 
construct the transformation T : Ti ^ L2 (cf. ( ^.ig| )): it 
is given by a combination 

\Tt(a/5)Td(ln6) 



T ■ 



Vb 



of a (complex) translation 
Tt{a/b) = expi{a/b)pH , {Tt{a/byj){s) 

and a (complex) dilatation 



■0(S- 



7d(ln6) = cxp y-{xnPH +PnXH)^^^b 
(Td(ln6)i/;)(s) = %/5?A(6s) . 



(2.40) 

a/b) , 
(2.41) 

(2.42) 
(2.43) 



III. DISCUSSION OF VARIOUS MODELS 

Unfortunately, most works dealing with non- Hermitian 
Hamiltonians are based on numerical methods which are 
not very well adapted for studies of the structure of the 
space Ti.. But note that an investigation in this direction 
has been started in We will here consider exactly 

solvable models. Among those works that also treat ex- 
actly solvable models, we mention [|l3| , where a similar 
construction of a Hilbert space as the one outlined in Sec- 
tion H has been recently applied to the study of certain 
coherent states. 

In the following we will throughout suppress the sub- 
script Ti for operators that act in the space Ti.; for oper- 
ators acting in L2 we keep the subscript i2- 



A. The one-dimensional Coulomb problem on the 
real half axis 



Let us start with a very simple example. We first dis- 
cuss it in a way inspired by Bender et al. and then 
reinterpret it according to the framework of Section ||. 

Consider the one-dimensional Coulomb problem on the 
real half axis a; > 0: The Hamiltonian reads 



Hl\a) = ipl?-^ 



(3.1) 



the eigenfunctions for real a; > and a > that are finite 
at the origin, have a finite slope there, and are decreasing 
for X —> 00 are ipotl 



1>„(^;a)-c-"-/(2")-Li^2i(ax/n) , 



(3.2) 



n = 1,..., with the Laguerre polynomials L^^-i- The 
eigenvalues are 



(- 

\2n 



(3.3) 



It is clear that $„(a;;a) is normalizable in ^2(0, 00) for 
real and positive a; the spectrum in this case is a bound- 
state spectrum. 

We may now use the solutions (|3.2| ) with a = ia, a real 
and positive, 

Mx) = Mx;ia) = e-'"^/(2")l^iW^(iaa;/n) . 

(3.4) 

The energies are then still real, 



- 2n/ 



and the Hamiltonian is 



TT 2 1" 

H = p 

X 



(3.5) 



(3.6) 



thus non-Hermitian with respect to L2. The asymptotic 
behavior of '(pnix) for complex x = |a::|e"' and large |x| is 
given by 



\i^n{x)\ ~ exp ( — — sini? 

' In 



(3.7) 



hence ipn converges most rapidly on the anti-Stokes line 
given by = 37r/2. The eigenfunctions (3.4) have no 
finite norm in ^2(0, 00). One would therefore consider 
the ipn{x) on the anti-Stokes line 



x{s) 



< s < 00 



(3.8) 



In terms of our general discussion in Section p B| , 
we would define the scalar product in the space TL by 

(cf. dii)) 



5 



H 



ds [V'(e"'s)]*V'(e'*s) 



(3.9) 



This corresponds to the special case (2.37) with a — 



and 6 = e . The canonical operators are given by (2.2t) 



and (2.28) 



e'V 



(3.10) 



therefore the canonical form of the Hamiltonian (3^) is 



\p^Y-io-^^ =-{fY+ (3.11) 



and is obviously Hermitian. Its L2 realization via T 
H ^ L2, ^{s) = (T^PKs) = V'(e'''s), is 



^L2 



and the transformed eigenfunctions read 
V3„(s) = $„(s;a) . 



(3.12) 



(3.13) 



In summary, the quantum-mechanical interpretation 
based on the canonical form (3.11) is completely equiv- 
alent to the original, Herm itia n problem; the change of 
the sign of the spectrum ( |3.5| ) is reflected by the sign 
in ( 3.12 ). Hence the solutions (|]J) do not contai n any 
new information. Moreover, the kinetic energy in (3.11) 
has — contrary to ( |3.6| ) — the wrong sign. This makes the 
model unphysical. 



B. The model of Cannata et al. 

In the model introduced by Cannata, Junker and Trost 
LSI the Hamiltonian is 



(3.14) 



(a similar model has been considered in |21|]). The Schro- 
dinger equation for this Hamiltonian can be reduced to 
Bessel's differential equation, so that its general solution 
is given by 



vE',(x)=cii7«(e'^) + C2i?(')(e'^ 



(3.15) 



For real x no normalizable solutions exist. But if x is 
chosen along curves in the complex plane, different pos- 
sibilities occur. We briefly summarize the main cases; 
details can be found in [[l8| . 

In the upper-half x plane, one can find normalizable 
solutions for every real value of the energy, i. e. no quan- 
tization condition appears and the spectrum has no lower 
bound so that the system is unphysical. If one divides 



the lower-half x plane into vertical strips separated by 
the lines Re x = mr with integer n, and if one assumes 
that the curves along which the wave functions are de- 
fined have two asymptotes that are vertical lines going 
down, then the following configurations occur: If a curve 
starts and ends in the same strip, no quantization con- 
dition appears, the spectrum is again unbounded from 
below. If the curve starts in one strip and ends in the 
next strip, then no normalizable solutions exist. If, how- 
ever, the curve starts in one strip and ends in the next- 
to-next strip (e. g. one asymptote of the curve is the anti- 
Stokes line at Rea; = 7r/2 and the other asymptote is at 
Rex — 57r/2), the quantization condition v = n + 1/2, n 
integer, appears. One then has the solution 



(n + 1/2)2 



(3.16) 



In addition to this case, on which our subsequent discus- 
sion concentrates, there exist other curves with similar 



spectra |18 



It follows from ( 3.16 ) that 

V'n - i(-l)-"-V-. 



(3.17) 



This means that the ipn with n > constitute a set of 
linearly independent eigenstates, so that one can start 
from the space V = spanl-00, i^i, ■ ■ ■} to construct H. 

The recursion relations zdZi,/dz ± vZ, 
Bessel functions can be used to find 



±zZu^i for 



1 



2n 
1 



1 



{■4'n 



n+1) 



(3.18) 
(3.19) 



Thus the two operators e^'^ and e~'^p are well-defined 
operators H ^ H. The two relations (3.18) and ( 3.19| ) 
lead to the version 



[{n + 1/2)2 _ ^ 



(3.20) 



of the Schrodinger equation Hipn — Enipn- (Note that 
these are algebraic relation s tha t hold for all Bessel func- 
tions and all n € C. It is ( 3.17 ) and the definition of V 
that are specific to the solution (3.16).) 

We first consider the question whether one can define a 
scalar product along a curve such that the Hamiltonian 
(3.14) becomes Hermitian. Since ( ^.14 ) is of the form 
(2.35), it is clear that only straight lines x = a + bx'^ 
are possible. Then p'^ = bp, and the Hamiltonian can be 
written as 



i^(x^p^) 



2&2 



(3.21) 



If one wants this to be Hermitian, only lines with Re a = 
mi/ 2 and Re6 = are allowed. Although such lines may 
be the asymptotes of curves that lead to the solution 
( |3.16 ), one cannot get a connected curve by joining these 
lines. 
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Ther efore , we now return to the general discussion of 
Section [I A and define the Hilbert space H by imposing 



the ort hogo naUty condition (ij jn, il' m.)H = for n,m > 
0, see (2.8). The equations ( 3.1g ) and ( 3.1£ ) can then 
be used to calculate the adjoints of e~"^ and e~"^p; the 
results are 



i>7i 



1 



and 



2(71 - 1) + 1 



1 



2(n+ 1) + 1 



P 



V'n+i (3.22) 



(3.23) 



It follows that p"^ ~ ie '^p is Hermitian. Furthermore, 
[x,p] = i leads to [e'^,p'^] = — i. Hence the two operators 



= ie-> (3.24) 



are Hermitian and fulfill = i so that they form a 

canonical set. 

But it seems that one cannot write the Hamiltonian 



(3.14) only as a function of i'^ and fp . One can readily 
express H as H — H{Sf ,p'^ ^x"^'), where x'^' — i(e'^ — 
e~'^ )/2 is another Hermitian operator. This operator 
fulfills [i'^'jp'^] = 0. This may mean that x'^' is only a 
function of ^ but we were not able to show this. 

One therefore has to find other canonical variables. We 
emphasize that for the following discussion we use the 
term Hermitian in the precise sense of symmetry^ i. e. an 
operator A is called Hermitian if Aip')-}i — {AiP,i/j')t-c 
for all ip, -0' in the domain of definition of A. 

We ca n use the method described at the end of Sec- 
tion [I A to try to find a new scalar product that makes H 
and e^'^ simultaneously Hermitian. As such a new scalar 
product is connected to the old one by (V'n,'0m)>^ = 
7n(V'n, V'm)-H, 7n > (scc ( 2.14 )), a Hermitian e~'^ would 
require real matrix elements {ipn, &~^^'ipn)H for all i^- But 
since ■(/'-i = iV'Oj one obtains ('(/'Oj e~'^V'o)'H = i; so that 
one cannot construct a scalar product for which c^'^ is 
Hermitian in the entire space Ti.. The best one can do 
is to make e~'^ Hermitian in the subspace spanned by 
{ipi, 7/'2, • • •}■ This is achieved by defining 



= v^^n = Vn + l72'/'«/v2 for n>0 , (3.25) 

and introducing a new Hilbert space Ti. with the new 
scalar product (V'n,V'ni)?^ ~ n,m>0. Then 



1 



1pn-l 



'n+1 



2v^n+ 1/2 ^v^n- 1 + 1/2 ^n+l + l/2j 

(3.26) 

for n > 1 and 



^^0 



. 7 , "01 



(3.27) 



These equations show 

(^„, e-'^V^™)^ = (e-'^V'n, (3.28) 
for all n, m > except for n = to = 0, i. e. the restriction 
e-'^l^^ : Vi (3.29) 

of e~'^ to the space Vi span{ipi, . . .} is Hermitian. If 
we define 



n+1 



y^n-l + 1/2 y^n + l + 1/2 



for n > 1 



(3.30) 



then the image of Vi under e is given by 
span{^i, . . .}. Consider the series 



JV 

■fi^^ = y^(~i)''^"+2fc+i 

fe=0 
tpn 



v/w + 1/2 
It fulfills 



+ (-1) 



v/[n + 2(7V + l)] + l/2 



H 



(3.31) 
(3.32) 

(3.33) 



^[n + 2(Ar + l)] + l/2 
and therefore converges in TL, i. e 



for TV ^ cx) 



V^^+T/ 



r7^=E(-l)'f-+2fc+i for">0 , (3.34) 



fc=0 



which means that the linearly independent (although 
non-orthogonal) vectors n> 1, lie dense in the space 
Ti. We now define the operator 



This operator acts as 



(3.35) 



x^in = 2Vn + l/2i^„ , (3.36) 
has a dense domain of definition in Ti and is Hermitian: 

{.^n-jX ^m^'fi ~ '^{^n.m+l ^f^n+l.rn) ^ {x ^njCm)^^ (3.37) 

for all n, m > 1. We use x'^ as one of our canonical 
variables. Its conjugate momentum can be defined by 



P 



-(e '> + pe~''') 



(3.38) 



The operator x'^ has to be distinguished from the op- 
erator e'^. Furthermore, although x'^ is well-defined in 
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the space by ( 3.36 ), in acting on the eigenfunctions 



tjjn via ( 3.34 ) it leads to the divergent series 



x'^ipn = 2^" + 1/2 

OO 

X ^(-1)''V" + 2/c + 1 + 1/2-0 



(3.39) 



fe=0 



so that x'^^lJn cannot be d efined by this series. Because of 
this, the construction of ( ^.38 ) and the demonstration of 



the commutation relation [ir^ p'^] = i is a rather techni- 
cal problem. The corresponding proofs are described in 
detail in Appendix Here we just note that 



2^n+ 1/2 -1 + 1/2 + 1 + 1/2^ 

(3.40) 

for all ri > 0, and that indeed 

{in,firn)i:i = {f£.n,irn)-H for aU n, m > 1 (3.41) 

and 

[i'.nfn^ien foraUn>l . (3.42) 

Now we are finally in the position to write down the 
canonical form of t he H amiltonian of the model of Can- 
nata et al.: From ( AlO|) we have p — iSfp^ + 1/2, and 
as explained in A ppen dix one can replace e'^ by x'^ in 
the Hamiltonian ( 3.14 ), so that H can be written as 



H = H{x-,f) = - 



-(x^f{p^f + 2ix^p^ + - + {x^f 



(3.43) 



One thing should be noted here: Although ( 3.14|) mul- 
tiplied with %l)n is of course well-defined in 7i, the two 
terms p^ipn and e^'^V'n are diverge nt; t heir divergencies 
cancel in the sum. With respect to ( 3.43 ) this means that 
one can group the terms on the right-hand side in such a 
way that no divergencies appear in intermediate results. 
This is achieved by writing (see (A13)) 



2 



x^ [x^ (1 - {ff ?) + 2if] + i 



(3.44) 



This form of the Hamiltonian is an operator that maps 
V = span{?/)o, "01, • • •} to itself. In this sense, ( 3.43| ) is 
equivalent to ( 3.14 ) in the entire space V. 

Let's consider the realization of this Hamiltonian in 
i2 (-00,00), i. e. = H{xl^,p'iJ, in which xl^ and 
p'j^^ are identified with the canonical L2 operators: 



2 „ d 1 

dx^ da; 4 



(3.45) 



The general solution of the Schrodinger equation Hl^ 
^ + l/2)2$^/2 is given by || <^>^{x) = Ci^^(x) + 
C2(p-u-i{x) with 



1 



ip^{x) = — J^+i/2(a;) 



(3.46) 



For the case corresponding to the spectrum ( 3.16 ), 
namely integer v — n, these solutions are proportional to 
spherical Bess el fun ct ions cij„(x ) + C 2j-n-i{x). In fact, 



the equations ( t3.26[ ), ( |3.27D and (|3.40 ) are closely related 
to the recursion relations for spherical Bessel functions 

The solutions (|1|) fulfill (2|] 



(3.47) 



for Re(i'* + /i) > —1. Here the value of ipy{~'x), a; > 0, 
has been chosen above the cut, i. e. —x = e^'^x. This 
result has quite interesting consequences: First of all, for 
integer n > only but not ip-n-i is L2 normalizable. 
Furthermore, the functions 



i(x) = ^Jn+ 1/2 (^„(x) = y ^"^ 3n{x) , 



are orthonormal in L2'- 



n > 0, 
(3.48) 

(3.49) 



in other words, the linear transformation T : 7i — > L2 
defined by Ttpn = f'n is an isometry. But it is not a 
unitary transformation, because the do not span the 
whole space L2. This can be seen most easily by consid- 
ering 



Dix,y) = 2^ ipn{x)ipl{y) 

n=0 



271 + 1 



n=0 



in{x)jn{y) 



(3.50) 



Completeness of the ipn would require D{x, y) = 6{x — y), 
but instead one obtains [GSl 



1 



D{x, x) — — 



(3.51) 



Even worse, the formally L2 Hermitian Hamiltonian 
(3.45) is not even Hermitian: It follows from (3.47) that 



^ — 7r In 



l|2 



sinh 27r Ini i 



Rei^+ 1/2 27rlmi/ 



for Kev > 



1 

'2 ' 
(3.52) 



thus H]^.-^ has eigenfunctions with complex eigenvalues 
E,, = (i/+ 1/2)2/2 e C in L2, so that {ip^,HL2^p^)L2 = 
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Our initial aim was to appl y the uniqueness theorem 
to the canonical formulation ( 3.43 ). But we have not 
checked whether the canonical operators Sf and p'^ are 
self-adjoint (as required by the uniqueness theorem). Al- 
though the Hermiticity and dense domain of definition of 
these operators are a necessary condition for this [0, it 
is in general not an easy task to check the self-adjointness. 
In our case, a simple argument can be given that shows 
that no unitary transformation T : 7i — s- L2 that maps 
x'^ and ff to the canonical position and momentum op- 



erators of the space L2, x'l 



Tx^T- 



exists: The operator p'^ is bounded, from (3.4C) one finds 



Tp^T- 



for all e H 



(3.53) 



Under a unitary transformation T the boundedness prop- 
erties do not change, hence the bounded operator jf can- 
not be mapped to the unbounded operator . 

This argument can even be generalized: The operator 



p" in ( 3.38D is not the only momentum operator that can 
be assigned to x'^, because if f{x'^) is a Hermitian function 
of a;'^, then 



P 



(3.54) 



is also a Hermitian operator with [i'^,7r'^] = i. Now sup- 
pose that this set {x'^ , tt*^) can be mapped to the canoni- 
cal position and momentum operators of L2 by a unitary 
transformation T. Then for ip £ L2, Eq. (p.54|) gives 



{pl,p){x) = (Tp^T-V)(x) = -i^ - fixMx) (3.55) 



da; 

= {Upl^U-'^){x) 



(3.56) 
(3.57) 



where dF/dx = f{x), and where the unitary operator U : 
1/2 L2 is the gauge transformation U = ey.Yi\F{x\^). 
Thus the canonical L2 momentum would be given by 



^u-'pI^u 



(3.58) 



Here all operators on the right-hand side are bounded, 
in contradistinction to the fact that is not bounded. 
This shows that if p'^ is bounded, then also all operators 
tt'^ of the form (3.54) cannot be unitarily equivalent to 



the canonical L2 momentum. 

Physically this means that one cannot interpret x"^ and 
p'^ or n'^ as position and momentum operators, resp. Thus 
our conclusion is that the physical interpretation of the 
model of Cannata et al. is still unclear. Of course, it may 
be possible that other, self-adjoint canonical operators 
can be found. As we have already pointed out, such 
a formulation is, in our opinion, necessary for the very 
definition of the model. 



C. The model of Hatano and Nelson 

In the context of studies of delocalization phenomena, 
the model of Hatano and Nelson Q has attracted a lot 
of interest recently ||2^,|2^ . It is defined in one dimension 
by the non-Hermitian Hamiltonian 



H{9) 



{P + 1.9)' 
2m 



Vix) 



(3.59) 



where 5 is a real parameter connected to an externally 
applied magnetic field, and V{x) is a random potential. 
It has been demonstrated numerically that at a certain 
critical value g = gc sl localized wave function turns into 
a delocalized one, and it has been suggested that this 
behaviour signals the occurrence of a delocalization phase 
transition. 

The numerical demonstration of a delocalization tran- 
sition at the critical value gc is based on the use of a mod- 
ified scalar product in the space of quantum-mechanical 
states . There has been some controversy about this 
point |2q| , but at present a consensus seems to have been 
reached that the delocalization transition is only visible 
if the density distribution of particles is calculated ac- 
cording to a scalar product based on the pro duct of left- 
and right-eigenfunctions of the Hamiltonian ( 3.59 ) p7|| . 

We want to point out here that a physical interpreta- 
tion that is based on this modified scalar product is only 
consistent if the underlying Hilbert space H is chosen ac- 
cordingly. In fact, as we will show, this Hilbert space is 
just of the kind c onsid ered in Section ||, i. e. below gc 
the Hamiltonian (3.59) is Hermitian in Ti.. But it turns 
out that under reasonable conditions this Hilbert space 
explicitly depends on the parameter g, so that the g de- 
pendence of the Hamiltonian ( |3.59| ) is not the complete 
dependence of the model of Hatano and Nelson on g. 

It is the aim of this Section to investigate this depen- 
dence. This is important for the phenomenological im- 
plications of the model of Hatano and Nelson, because 
g is considered to be an external parameter that can be 
varied. 

We want to emphasize that we are dealing here only 
with the quantum-mechanical model of [0, not with its 
applications to statistical-mechanical problems like vor- 
tex depinning in type-II superconductors or population 
biology (2|. 

Consider a set of eigenfmictions 'ipn{x;g) of ( 3.59| ), 



H{g)ipn{x;g) = En{g)i)n{x; g) 



(3.60) 



determined as solutions of the corresponding differential 
equation and being elements of a Hilbert space Tip (e. g. 
the space L2{—oo, 00) or the space of periodic functions). 
In , examples are given in which the eigenvalues En (g) 
are real for g being smaller than some critical value gc, 
for g > gc they become complex. 

Let us first recall the construction of the modified 
scalar product as introduced by Hatano and Nelson. We 
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refer here to the detailed description given in |27 
cording to which the ^/;„ can be chosen such that 



dxiJnix;-g)*ipmix;g) = [^n{-g),'4'm{g))na (3-61) 



is always well-defined and given by 



(V'«(-g),'0m(ff)) 



Ha 



(3.62) 



The probabilistic interpretation of the model is based on 
this bilinear form; the integrand in ( 3.6l| ) for n = m, 



q{x) = 'ipn{x;-g)*^n{x;g) = .g) (in the 

notation of |2^, the superscripts L and R denoting left- 
and right-eigenfunctions) is interpreted as the particle 

density in the bulk of the sample. 

It is important to realize that ( |3.62 ) rests upon some 
quite non-trivial assumptions. First of all, as t he Tip 
adjoint of H{g) is given by {g) = H{-g), Eq. ( |I6^ ) 
can only hold provided the spectra of H{g) and H{—g) 
can be ordered in such a way that [p7| 



E^i^g)* = E^{g) 



(3.63) 



By taking the complex conjugate of the differential equa- 



tion corresponding to the Schrodinger equation (3.60) 
we see that if H{g)ipn{x;g) = En{g)ipn(x] g) , then 



H{gWn{x\g) = E^{g)'iJj*Xx; g) 



(3.64) 



i. e. th e eig envalues come in complex conjugate pairs. 
Hence ( ^.63 ) asserts that the spectra of H{g) and H{—g) 
coincide. 

In order to construct a new Hilbert spa ce th at is 
based on a scalar product consistent with ( 3.62 ), we 
again start with the infinite-dimensional vector space 
V{g) = span{'i/',x(5), ft = 0, 1, . . .} of the finite superpo- 
sitions of the eigenfunctions. Next, we define a linear 
operator Mg : V{g) V{—g) which acts as 



Mgil^nig) = ipn{-g) 



(3.65) 



and use it to define a scalar product in V{g) through 



■Ho 



(3.66) 



Proceeding as in Section O. 
with respect to the norm^lV'll v(g) 



we can complete this space 
\/(V'7'0)v/ol" to ob- 



tain the separable Hilbert space Ti.. Then, by (3.62), the 
eigenfunctions ipnig) form an orthonormal basis in Ti.. 

There may also be cases in which one is co nsidering 
an enumerable set of solutions ipnix'^g) of (3.60) that are 
not elements of a space Ho (e. g- because they are not 
normalizable in Tio)- One can then still construct TL in 
complete analogy to Section || by considering this set of 
solutions ^n{x',g), defining V{g) as above, and defining 
the scalar product in V{g) so as to fulfill 



('0n(g),'0m(.g))-H = 5n 



(3.67) 



This definition is always possible; it does n ot re ly on 
additional assumptions like (|3.63 ). Condition ( 3.63|) only 
has to do with the explicit form of the metric operator 
Mg relating the scalar products in Ti. and Tin . 

We can now formulate the model of Hatano and Nelson 
entirely in the space Ti, and since the spectrum is real 
for g < gc, we find as in (2.10) 



HHg)^H{g) for g<g. 



(3.68) 



We emphasize that it is essential for the following dis- 
cussion that the metric operator Mg in (3.65) depends 
on g, because then it may happen that the space Ti will 
be different for different values of g. 

For matrix elements of an operator A-^ and vectors ip, 
■0' that are contained both in Ti. and Tip, and for which 
both the Hermitian adjoint A\_^ in Tio and in Ti are 
defined, one finds 

(A*,7/-,V')h = ii',Ani^')n = (MgV,A„V')Ho 
= iAl^Mgi;,,p')n„ = (M^M-g^j^Mgi^, V')ho 
^ iM^gAl,MgiJ;,^')H (3.69) 



since MgM-g = 1, so that 



Al^ = M.^AIm, 



(3.70) 



Although Ti is well-defined by the construction out- 
lined above, we do not know much about it as the eigen- 
functions ipn{x; g) are not explicitly known. (To our 
knowl edge, no non-trivial, exactly solvable model of the 
form ( 3.59|) that exhibits a delocalization transition ex- 
ists.) One therefore has to rely on some assu mptio ns, 
and we will assume here that the operator x in ( 3.59| ) is 
Ti Hermitian, 

x^^x . (3.71) 

This assumption is quite natural, because according to |^ 
the fact that a delocalization transition (leading from a 
real spectrum to a complex one) occurs, does not depend 
on the detailed form of the potential V{x). Since below 
gc the Hamiltonian H{g) is Hermitian (cf. ( |3.68| )), and 
since this has to be true for a large class of potentials, the 
operator x should be Hermitian. Note that here , cont rary 
to the treatment of the examples in Sections HI A and 



HIB, we treat a; as a physical observable. 

Consider the case g < g^- We then have from the Her- 
miticity of H{g) and x 



{P^ - i.g)' = (P + i.9)' 
This equation has two solutions: 

■^^ — and 



p' 



-P 
p + 2ig 



(3.72) 

(3.73) 
(3.74) 



The first solution ( 3.73 ) would be compatible with an 
Mg that is independent of g, sin ce M^gpMg = —p in this 
case. The second solution (3.74), 
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M_gpMg =p + 2ig 



(3.75) 



requires a g dependent Mg. But the first sohition is 
physically unacceptable. The reason is that (as H{g) 
is Hermitian) we can assume that the time evolution is 
governed by the Schrodinger equation, so that the expec- 
tation value {x){t) = {'ip{t),x'i{:{t))-H evolves according to 



-{x)it)=i{[H{g),x]){t)^{v)it) 



(3.76) 



with V = (jj + ig)/m. Hence p -\-\g = p is the kinetic 
momentum. (This observable has also been considered in 
|0|.) The first solution (3.73) then corresponds to an anti- 
Hermitian velocity operator, = — w, giving a purely 
imaginary expec tation value for the velocity. Only the 



second solution (3.74) leads to a real expectati on va lue 
for the velocity. Therefore, we only consider ( 3.74 ) to 
be physically meaningful, and consequently Mg has to 
depend on g. 

We now want to discuss how the quantum-mechanical 
predictions of the model vary with g. The situation in 
quantum mechanics is similar to the situation in classical 
Hamiltonian mechanics: There one has two canonically 
conjugate variables x'^ and p^ which are the independent 
variables, and the Hamiltonian has to be expressed in 
terms of these two variables: 



H = H{g-x^,p-) 



(3.77) 



Varying any parameter g appearing in the Hamiltonian, 
the system reacts in precisely the w ay given by the g 
dependence of the formulation (3.77). In the quantum- 
mechanical case, these two canonically conjugate vari- 
ables x"^ and p"^ have to be two operators fulfilling canon- 
ical commutation relations [a:;'^,p'^] = i and have to be 
self-adjoint. By this last requirement, the structure of 
the underlying Hilbert space comes into play. Having 
such a formulation, we can use the uniqueness theorem to 
map H{g;x^jp'^) uniquely (modulo unitary equivalence) 
to the L2 formulation Hl^ = H{g; x'i^,p'i^). Thus, the 
g dependence of H{g; x^^^pl^) describes the complete g 
dependence of the physical system. 

Applying this reasoning to the model of Hatano and 
Nelson, one has to find the two canonically conjugate 
operators x'^ and p'^. Again, we assume that x is Hermi- 
tian in H, whereas p"^ is non- Hermitian, see (3.74). But 
as [x,p\ = i, one finds that one may choose 



x'^ — X and p'^ = -^{p + P^) 



(3.78) 



as the set of Hermitian operators wi th ca nonic al com mu- 
tation relations. (For g < g^ Eq. ( [3.7g ) and (3.74) just 
give p"^ = p + ig = p'^'" . ) 



The next task is then to express the Hamiltonian (3.5£ ) 
through x'^ and p'^. In order to do this, one has to know 
how p^ acts in Ti. for arbitrary g. Again, this would be 
clear once we knew the states ipn explicitly. But since we 



don't know them, we have to rely on some assumptions, 
and here we will restrict our discussion to a physically 
rather uninteresting case, namely the case of wave func- 
tions that are defined with infinite, non-periodic bound- 
ary conditions. It is well-known |^ that under these con- 
ditions there is no g dependence at all. The numerical 
demonstrations of delocalization transitions have always 
been carried out for wave functions with periodic bound- 
ary conditions. Still, the example of infinite, non-periodic 
boundary conditions serves illustrative purposes. At the 
end, we will comment on the case of periodic boundary 
conditions. 

The point is that for non-periodic boundary conditions 
the unitary transformation Tg : Ti. L2 is known to be 
only a function of x ~ x'^ 0]: 



(3.79) 



Such a Tg acts on states tp ^ H simply as {Tg{x'^)ip){x) 
Tg{x)il](x), so that 

(^,^')w = {Tgi,,Tgi^')L, = / dx\Tg{x)\''r{xW{x) 



(3.80) 

From this equation one can deduce how p^ acts on 
states of Ti: A simple calculation based on (p^ip^ 4'')'H — 
{'4',P'>P')n and partial integration shows that 



(3.81) 



Thus, from (|3.78| ) we obtain 



(3.82) 



Substituting this into ( 3.5S| ) leads to the expression 

H{g)^H{g;x'=,p') 
W . i d 



Vix") . (3.83) 



Equation (3.83) represents the canonical form of the 
Hamiltonian of the model of Hatano and Nelson for in- 
finite, non-periodic boundary conditions. As before, in 
this form the model can immediately be formulated in 
the space L2 by just replacing x'^, p^ by x\^, p^^, resp. 
As one can see, we have canonically coupled a nontrivial, 
purely imaginary field 



A{x;g)^--\n\Tgix)\'+,g 
to the system with 5 = 



(3.84) 



The spectrum of ( 3.83 ) is real for g < gc and complex 
for g > gc- This means that H{g;x^,p'^) i s He rmitian 
below gc and non-Hermitian above gc (see ( |3.68 )). The 
Hermiticity translates into a vanishing anti-commutator 
{p'^,A{x'^;g)} = for g < gc- Formulating this in L2, 
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{-id/dx, A{x; g)}(p{x) = for all (p € L2 , (3.85) 

leads to the equation Lp{dA/dx) = —2A{dLp/dx) which — 
as it has to hold for all ip G L2 — requires 



A{x;g)=0 (for 5 < gc) 



(3.86) 



For g > gc the spectrum becomes complex, and A{x; g) 7^ 
is possible. Hence one can write explicitly 

H{g; x^,p^) = -^b^ + 6(5 - gc)A{x^;g)? + V{x^) . 
Zm 

(3.87) 

In fact, for infinite, non-periodic boundary conditions 
it is well-known [Q that Tg is given by 



leading to 



Tgix'') = exp(-5a;'' 



A{x; g) = for all g 



(3.; 



(3.89) 



i. e. gc = 00 in (3.87), and the spectrum is real for all g. 

The transformation Tg in ( 3.88| ) is often refered to as an 
"imaginary gauge transformation" . Our analysis shows 
that this notion is quite misleading: Nothing has been 
gauged away by Tg] we just have two formulations of 
the same physics, realized in the two different, though 
unitarily equivalent Hilbert spaces Ti. and L2. 

The problem in the case of Hilbert spaces Hq and Ti 
that contain wave functions with periodic boundary con- 
ditions is that in these spaces only the differentiation 
operator p is defined; the operator x has a vanishing do- 
main of definition. This makes it impossible to construct 
a canonical formulation for which the uniqueness theo- 
rem m ay be used. Still, if one chooses a peri odic p oten- 
tial in ( 3.59 ), our reaso ning a pplies up to Eq. ( 3.75 ); only 
the argument given in (3.76) can — strictly speaking — no 
longer be used. This means that also in the case of peri- 
odic boundary conditions it may happen that the space 
Ti. is an unknown function of g, making it a priori impos- 
sible to compare the predictions of the model for different 
values of g. 

To our knowledge, this aspect has not been discussed 
in the literature up to now, but it is in our opinion of 
fundamental importance for the phenomenological inter- 
pretation of the model of Hatano and Nel son. W e see, for 
example, that the canonical formulation ( 3.87 ) describes 
a non- continuously differ entiahle change of the Hamil- 
tonian at the transition point gc. Whether such a be- 
haviour should be called a phase transition is not clear 
to us; it seems to be more similar to the sudden switching- 
on of a magnetic field as the parameter g passes gc- It is 
clear that the behaviour of the system will then change, 
localized states may become delocalized, but one would 
hardly call this a phase transition. 

Let us mention another aspe ct of the model of Hatano 
and Nelson. We have seen in (3.6S) that as long as the 



spectrum of H{g) is real, the non-Hermiticity of H(g) 
with respect to the space Tio is a rather superficial one 
that can be avoided by considering the space H. Thus 
the time evolution can consistently be described by the 
Schrodinger equation, leading to a unitary time-evolution 
operator. The situation changes when the spectrum be- 
comes complex. It is then no longer possible to con- 
struct a Hilbert space in which the Hamiltonian is Her- 
mitian and which at the same time contains eigenstates 
with complex eigenvalues. Although a completely con- 
sistent probabilistic interpretation is still possible, the 
interpretation of the Hamiltonian as the generator of the 
time evolution would then lead to a non-unitary time- 
evolution operator with all the problems mentioned in 
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IV. CONCLUSIONS 

In this work we have studied the quantum-mechanical 
interpretation of models with non-Hermitian Hamiltoni- 
ans (in the usual sense of the space L2) and real spectra. 

Assuming that the systems under consideration are 
closed, and that the Hamiltonians are the generators of 
the time evolution for these systems, we construct sepa- 
rable Hilbert spaces in which the Hamiltonians are Her- 
mitian. 

Within this construction, we set up a canonical for- 
mulation in which the Hamiltonian is only a function of 
two Hermitian operators x'^ and p'^ that fulfill canonical 
commutation relations. Then a unique physical inter- 
pretation is obtained if these operators are self-adjoint, 
because in this case the model can be formulated as a 
Hermitian problem in the space L2. If, on the other 
hand, this equivalent L2 description cannot be achieved, 
we consider the model to be quantum-mcchanically in- 
consistent. 

We apply the above construction to a number of mod- 
els recently discussed in the literature. 

As a first example, we have analysed a simple non- 
Hermitian model that is based on the Hermitian one- 
dimensional Coulomb problem on the real half axis. Mo- 
tivated by the approach of Bender et al. we extend 
the eigenfunctions to complex values of the coordinate. 
We then show that the canonical formulation of this ex- 
ample does not contain any new information compared 
to the original, Hermitian version of the model. 

The next example is the model introduced by Cannata 
et al. [Q. Although here we also find a canonical for- 
mulation, it turns out that no equivalent L2 formulation 
exists. In our opinion this may indicate that the model 
is intrinsically inconsistent. 

Our last example is the phenomenologically important 
model of Hatano and Nelson Here the appropriate 
Hilbert space is already determined by the introduction 
of the modified scalar product in j|] . Under these condi- 
tions, the dependence of the Hamiltonian on the external 
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parameter g that is responsible for the non-Hermiticity 
of the Hamiltonian does not fully describe how the sys- 
tem varies with g. For the special case of wave functions 
with infinite, non-periodic boundary conditions, we give 
the canonical L2 formulation which clearly exhibits this 
complete g dependence. It turns out that at the critical 
value g — gc, where the states of the system undergo 
a delocalization transition, also the L2 version of the 
Hamiltonian changes in a non-continuously differentiable 
way from Hermitian to non-Hermitian. In our opinion, 
this makes the usual interpretation of this transition as 
a phase transition questionable. We argue that in the 
phenomenologically relevant case with periodic bound- 
ary conditions similar effects may also occur. 
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APPENDIX A: PROPERTIES OF THE 
CANONICAL VARIABLES FOR THE MODEL OF 
CANNATA ET AL. 

In ( |3.35| ) the operator x'^ is defined on the subset 
{^1,^2, ■ ■ •} which is dense in Ti. Eq. (3.37) shows that 



x'^ is Hermitian. We can therefore use the fact that every 
densely defined, Hermitian operator has a closed exten- 
sion [^ , to extend the domain of definition of x'^ to some 
larger space than Vj. The problem is that x'^ipn is not 
defined by this procedure (cf. ( |3.39| )). We will therefore 
in the following define every operator on the space Vj, 
check its algebraic properties there, and will then implic- 
itly assume that we consider closed extensions where this 
is possible. 



To start, we re write (3.26) and (3.27) with the help of 
(1^3^ and (U|) in the form 



n+l 



2(71-1 + 1/2) 2(n-hl-hl/2) 



(Al) 



for n > 2, and e^'^Ci = 6/5 - E^i ("O'^C/c- H follows 
that although 



e "'x'^Cn = 6 for all n > 1 



one obtains 



i'^e '"^^n = ^„ only for n>2 



(A2) 



(A3) 



whereas i'^e"'^^! is divergent. 

In order to construct the canonically conjugate mo- 
mentum operator ff, we first note that ( 3.19| ) leads to 



n+l 



- 1 + 1/2 + 1 + 1/2/ 

(A4) 



for n > 1 and e '^pi/^o 



use pe — c '^p = — e 



iV)o/2 - i/)i/(2%/3). Then we 
together with (|A2D to get 
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p = e ^^px'^ — 1. Combining this with (A4) gives for all 
n > 1 



Kn = - 1 + 1/2)C„ - {2n + 1) 



n+1 



1 + 1/2 



(A5) 



This equation shows that p is anti-Hermitian: 
for all n, m > 1. Now define 



p^ = --(e-> + pe--) = -ie-'-(p-l/2) . (A7) 



This operator fulfills 

^ ^" " 2 n - 1 + 1/2^""' + n + 1 + 1/2^"+' 



(A8) 



for all n > 1. In the matrix element iS,mP'^£,m)-fy 
— i(Cn, e~''^(p — l/2)^m)^ (with n,m>l) the vector 



{p - l/2)e™ = (™ - l)Cm - {2m + 1) 



m+l 



y/m+l + 1/2 

(A9) 



never has a component proportional to ipo (see ( 3.30 )), 
because for m — 1 the first term vanishes. Then, as a 



consequence of ( |3.28[) , the operator e~"^ in {^n,'P''Cm)-fy 
can be shifted to the left side for all values n,m > 1. 
This, together with the anti-Hermiticity of p shows that 
p'^ is Hermitian in the entire space V^. 

The commutator [x'^ , p'^] can be calculated in a similar 
way: One has x'^p'^^^n — —ix''e~^^{p — l/2)^„. Si nce ( p — 
l/2)^„ never has a component, one can apply (A3), so 
that 



Applying ^„ to it, one should note that, for example, 
(i'^)^Cra = 2\/n + l/2x'^V'„ is no t defined. But the terms 
on the right-hand side of (A12) can be arranged in the 
following way: 



H 



[x' (1 - (p^)2) + 2ip= 



(A13) 



Then all intermediate results are well-defined in H. One 
obtains 



Hix',f)^n 

(n- 1-^1/2)3/2 



(A14) 



(n-l- 1 -t- 1/2)3A 



+ 1 



for n > 2. We emphasize that this result cannot be used 

(n) 

to check the eigenvalue equation for — limAr^oo fj^ 
(cf. (|3.31| )): The Hamiltonian is an unbounded operator, 
hence one cannot expect that 



H hm fP 



hm HfP 



(A15) 



In fact, limN^Qo H fj^^^ is divergent. But one can apply 



tpn directly to (A13). Again, all intermediate results are 
well-defined, for example 



x^l - iff) +2if] ^„ = 



n(n -|- 1) 
2^n + l/2 



for n > 1 , 
(A16) 



and one obtains the expected result H'ipn — {n + 
1/2)2 -0„/2. gygjj turns out tha t thi s result holds for 
= 0, so that the representation (A13) is correct in the 
entire space V — span{-0Oi "^ii • ■ •}■ 



APPENDIX B: THE CONNECTION BETWEEN 
REAL SPECTRUM AND VT SYMMETRY 



x'^p'^^n = -i(p - l/2)Cn for aU 71 > 1 



(AlO) 



Writing p'^ = — i(p-|- l/2)e^'^, one can use ( [A^ ) to derive 
p'^x'^^n = ^i(-P + l/2)Cn for n>l. Combining these 
relations leads to the desired result: 



[i^p'^]^,, = iC„ for all n > 1 



(All) 



In summary, the operators Sf and p'^ are Hermitian, 
densely defined operators in 7i that fulfill canonical com- 
mutation relations. 



One can now express the Hamiltonian ( 3.14 ) as a func- 
tion of Sf and ff: If we restrict the discussion for the 
moment to the subspace span{^2, Csj • ■ -li can replace 
e^'^ in (pl^) by {x^f. Then (|XlO|) gives 



2 



,c\2 



(A12) 



In [g. Bender et al. attribute the reality of the spec- 
trum of non-Hermitian Hamiltonians H to the VT in- 
variance oi H. In fact, many examples of non-Hermitian, 
VT symmetric Hamiltonians with real spectra have sub- 
sequently been found Still, also examples of non- 
Hermitian Hamiltonians with real spectra that are not 
VT invariant are known, see, e. g., |J|,|l|. A proof that 
VT invariance is related to the reality of the spectrum 
does not exist. Here we want to present some arguments 
against such a relation. 

Let us consider general properties of non-Hcrmitian, 
VT symmetric Hamiltonians. We put O = VT in the 
following and assume 



[H, 6] = 



(Bl) 



Since O is anti- unitary, we find for an H eigenstate V'-E„; 
HipE^ — EnipE^ {En being complex in general). 
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HQ^E^ = QH^E„ = QE^^E^ = E*Qi,E„ , (B2) 



showing firstly that QipE,^ is an eigenstate of H with 
eigenvalue E* and secondly that the eigenvalues of a O 
symmetric, non-Hermitian Hamiltonian always occur in 
complex-conjugate pairs. Thus, if the spectrum of H is 
non-degenerate, one has 

0V'B„ = const V'_E* ■ (B3) 

We want to emphasize that the Q symmetry of a non- 
Hermitian Hamiltonian H is not sufficient to ensure the 
reality of the spectrum of iJ . Only if one makes the 
additional assumption that ipEn is a simultaneous eigen- 
state of 8, i. e. QipE„ = 'ypEn, one can conclude that the 
spectrum of H is real. 

But it is important to note that the last assumption is 
crucial for the reality of the spectrum, and that there is 
no reason to believe that this assumption holds in gen- 
eral. To see this, recall the usual argument that is used 
to show that a Hamiltonian H and a further linear oper- 
ator A can be simultaneously diagonalized if [H,A] — 0: 
If Hipn — Enipm one has 

HAi^r. = Ai?Vn = AEn^n = Er^A^n , (B4) 

showing that Aip„ is an eigenstate with the same eigen- 
value as ipn- If now the spectrum of H is non-degenerate, 
one has Aipn = const ip„ , and if the spectrum of H is de- 
generate, one may still diagonalize A in the eigenspace 
corresponding to En, so that the eigenstates of H with 
eigenvalue En may still be superimposed to give eigen- 
states of A. If, however, A is an anti-linear operator. 



then (B4) gives HAipn = EnAipn, so that A-0„ is not an 
eigenstate with the eigenvalue Hence, in general one 
cannot simultaneously diagonalize H and an anti-linear 
operator A, even if their commutator vanishes. 
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